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We consider the one-chain term in a skeleton expansion for Euclidean QCD obscrvablcs. Focusing on 
the particular example of the Adler D function, we show that although there is a Landau pole in the 
coupling at Q 2 = A 2 which renders fixed-order perturbative results infinite, the Landau pole is absent 
in the all-orders one-chain result. In this approximation one has finiteness and continuity at Q 2 = A 2 , 
and a smooth freezing as Q 2 — > 0. 



In this talk I want to describe some recent 
work with Paul Brooks ^ in which we con- 
sider the low-energy behaviour of Euclidean 
QCD observables. We investigate the Q 2 - 
dependence of all-orders perturbative resum- 
mations obtained using the so-called leading- 
b approximation l ^ l vrM where b — (33 — 
2Nf)/6 is the leading beta- function coeffi- 
cient in SU(3) QCD with Nf active quark 
flavours. This is closely related to the one- 
chain term in a skeleton expansion in which 
a single chain of fermion bubbles is inserted 
into a basic skeleton diagram We shall 
demonstrate that in this approximation the 
Landau pole in the coupling at Q 2 = A 2 is 
absent, and one has finiteness and continu- 
ity at this energy with a smooth freezing be- 
haviour at lower energies. We focus in this 
talk on the Adler Z?-function, but DIS sum 
rules are also studied in Ref. 
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The Adler D function is directly related to 
the vacuum polarization function W J,U (Q 2 ), 
(Q 2 = -q 2 > 0), 

n^(Q 2 ) = 167r 2 i / d?xe*- x (0\T(J„(x)J v (p))\0). 



Adler Function D(Q 2 ) 

d{Q 2 ) = -- A Q 2 ^(Q 2 )- (3) 

This may be written as a sum of the parton 
model result and QCD corrections T>(Q 2 ) 

#(Q 2 ) = 3]TQ 2 [1 + 2?(Q 2 )]. (4) 
/ 

These QCD corrections are split into a per- 
turbative and non-perturbative part 

V(Q 2 ) =V PT {Q 2 )+V NP (Q 2 ). (5) 

The PT component has the form 

??PT(Q 2 ) = a(Q 2 ) + $>na" +1 (Q 2 ). (6) 

n>0 

Throughout the talk we will take a(Q 2 ) = 
a s (Q 2 )/ir to be the one-loop form of the QCD 
coupling 

a{Q2) = 61n(Q 2 /A 2 ) • (7) 

The NP component in Eq.(5) will have the 
form 



Vnp{Q 2 )=Y, C \q2 



A 2 



(8) 



(1) 

Current conservation dictates that this has 
the tensor structure 

n^(Q 2 ) = ( q ^- gtw q 2 )n(Q 2 ). (2) 

Only n((5 2 ) — 11(0) is observable, so it is use- 
ful to eliminate the constant and define the 



The leading OPE contribution for the Adler 
function is the dimension 4 gluon condensate 
contribution 

1 

Q A 

We are interested in the behaviour of 
V{Q 2 ) = V PT {Q 2 ) + V NP {Q 2 ) as Q 2 -> 0. 



Go(a(Q 2 )) = ^(0\GG\0)C GG (a(Q 2 )) . (9) 
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Clearly at any fixed order perturbation the- 
ory breaks down at Q 2 = A 2 , the Landau pole 
in the coupling, and a(Q 2 ) — > oo. Evidently 
we need a resummation of perturbation the- 
ory to all-orders to address the freezing ques- 
tion, and we need to combine the resumma- 
tion with the OPE condensates. The large- 
Nf limit provides a way of formulating this 
resummation. The coefficient d n may be ex- 
panded in powers of Nt the number of quark 
flavours 



behaviour of the form d n ~ Kn 7 (b/2) n nl 
involving the QCD beta-function coefficient 
b = (33 — 2 AT/) /6, it is then natural to replace 
Nf by (33/ 2 — 36 ) to obtain an expansion in 
powers of &GEE] 

d n = 4") b n + dfr-Vb* 1 - 1 + . . . + 4°) (13) 



(L) 



4 n V 



d n = d^Nf + d^Nj- 1 



The leading large- Nf coefficient dn may be 
evaluated to all-orders since it derives from 
a restricted set of diagrams obtained by in- 
serting a chain of fermion bubbles inside the 
quark loop. The crucial ingredient in the cal- 
culation is the chain of n renormalised bub- 
bles Bfo, 



The leading-6 term d[ 
(— 3) n d n l ^b n can then be used to approximate 
d n to all-orders, and an all-orders resum- 
mation of these terms performed to obtain 
^[0] /jq\ 2?y^(Q 2 ). If we use the Borel method to de- 
fine the all-orders perturbative result we ob- 
tain 



(n) 



(k 2 g^ - fc M fc„) 
(fc 2 ) 2 



3 
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Here k is the momentum flowing through the 
chain of fermion bubbles. The factor in the 
square bracket is the one-loop vacuum polar- 
ization contribution Ilo(fc 2 ). The constant C 
depends on the subtraction procedure used to 
renormalise the bubble. With MS subtrac- 
tion C = — | . We shall choose to work in the 
"V-scheme" which corresponds to MS with 



the scale choice /x 



2 _ -5/3 



Q 2 , i n w hich case 



C = 0. The result for d [ n ] {V) is™ 



d [ :\V) = —{n- 



-2n 



2 n+2 



16 



§ + l>m>0 



m(l - 2' 



(1 - 2 2m -"- 2 )C 2m+ i] 



(12) 



This large- Nf result can describe QED vac- 
uum polarization, but for QCD the correc- 
tions to the gluon propagator involve gluon 
and ghost loops, and are gauge-dependent. 
The result for IIo(fc 2 ) is proportional to 
—Nf/3 which is the first QED beta- function 
coefficient, 6. In QCD one expects large-order 



dze- z/a{Q2) B[V ( p^](z)(U) 



The Borel transform B[V^]{z) can be found 
in RefPl It has the form 

B[v$}(z) = jr Mn) - Mn)z " 



A 1 (n)z n 



1 



OO 

£ 

n=l 



B (ri) +B 1 (n)z n Bi(n)z n 



0-t) 2 



(15) 



The expressions for the A 01 (n) and B 01 (n) 
are given m RefP. For the Adler function in 
leading-6 approximation we see that there are 
single and double poles in B[Dp^](z) at po- 
sitions z = z n and z = — z n with z n = 2n/6 
n = 1,2,3,.... The singularities on the pos- 
itive real semi-axis are the infrared renor- 
malons, IR n and those on the negative real 
semi-axis are ultraviolet renormalons, UV n . 
We shall see that they correspond to integra- 
tion over the bubble-chain momentum k 2 in 
the regions k 2 < Q 2 and k 2 > Q 2 , respec- 
tively. The IR renormalon singularities lie on 
the Borel integration path along the positive 
real axis. This leads to an ambiguous imag- 
inary part which is structurally the same as 
a term in the OPE expansion in Eq. (8) . The 
IR n renormalon ambiguity is in one-to-one 
correspondence with non-logarithmic UV di- 
vergences present in the (A 2 /Q 2 ) n term of the 
OPE SI OPE ambiguities and perturbative 
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ambiguities can cancel once a definite regula- 
tion of the Borel integral , for instance prin- 
cipal value (PV), has been chosen, and the 
PT and NP components are then separately 
well-defined. The (PV regulated) Borel inte- 
gral may be evaluated in terms of exponential 
integral Ei functions El but notice that the 
Borel integral diverges for Q 2 < A 2 , and po- 
tentially for Q 2 = A 2 ! 



Vf, e ^/«(Q 2 ) E i ( Zn \ 

h: Uq 2 )) 



(A (n) - z n A x (ri)) - z n A x (n) 



a(Q 2 ) 
(A (n) - z n Ai{n))\ 



E 

n=l 



a(Q 2 ) 



{B (n) + z„B x (n)) - z n B x (n) 

(16) 



a(Q 2 ) 
- {B (n) + z n B x (n))} 



This expression has the property that it is fi- 
nite and continuous at Q 2 = A 2 and freezes 
smoothly to a freezing limit of X>^(0) = 0. 
Similar behaviour is found for GLS / polarized 
Bjorken and unpolarized Bjorken DIS sum 
rules, K. ( p}(Q 2 ) and^(Q 2 ). SeeEfor their 
definitions. Interesting connections between 
these sum rules have been explored in RefP. 
The finiteness and continuity are delicate. 
The Ei functions potentially have a diver- 
gence proportional to lna(Q 2 ) as Q 2 — > A 2 , 
but the coefficient of this divergent term is 



z 2 n [A x (n) + B x (n)]. 



(17) 



rt+l 



For 1C^pt{Q 2 ) and U p L ^ (Q 2 ) the equivalent co- 
efficients are (-8 + 2 + 16-10 = 0) and 
(8 — 6 — 2) = 0, respectively. There is a rela- 
tion between IR and UV renormalon residues 
which ensures the divergent term vanishes 

z 2 n+3 B x (n + 3) = -z 2 n A x (n). (18) 



This ensures that 

oo 

J24lMn) + B 1 (n)] = 0. (19) 



71=1 



Another similar relation is El 

A (n) = -B (n + 2). 



(20) 



We shall show that these relations are un- 
derwritten by continuity of the characteristic 
function in the skeleton expansion. 

The one chain term in the QCD skeleton 
expansion can be written in the formal 



vf T {Q 2 



dt uj v {i)a{tQ 2 ) . (21) 



Here wji(t) is the characteristic junction. It 
satisfies the normalization condition 



dt u>t> (0 = 1- 



(22) 



t = k 2 /Q 2 , and so one is integrating over the 
momentum flowing through the chain of bub- 
bles. u>T>(t) and its first three derivatives are 
continuous ^ at t — 1, and the integral di- 
vides into an IR and UV part, corresponding 
to k 2 < Q 2 and k 2 > Q 2 , respectively. 



V 



(L) 



1 p OO 

rr , dtuj^ R (t)a(tQ 2 )+ / dtu% v a{tQ 2 ). 

(23) 

The first term involving uij^ reproduces the 
IR renormalon contributions, and the second 
term involving uj^ v the UV renormalon con- 
tributions. For Q 2 > A 2 the first integral en- 
counters the Landau pole in the coupling in 
the region of integration and requires regula- 
tion. If one uses a P V definition one can show 
by a change of variable that the one-chain 
skeleton expansion result is exactly equiva- 
lent to the PV regulated Borel integral of 
Eq.(14), and yields V [ p^(Q 2 ) as in Eq.(16). 
For Q 2 < A 2 the first integral no longer re- 
quires regulation, but the second UV term 
docs. If one uses PV regulation one can show 
that the one-chain result is exactly equivalent 
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to a PV-regulated modified Borel representa- 



tion 



1 10 



dze- z/a(Q2) B[V ( p}}{zl24) 



Note that the contour of integration now 
runs along the negative real axis in the Borel 
plane, and therefore there is now an ambigu- 
ous imaginary part due to the UV renor- 
malons. The ambiguity contributed by UV n 
is of the form (Q 2 /A 2 )™. These ambigu- 
ities are associated with IR divergences of 
dimension-six four-fermion ope rators associ- 
ated with UV renormalons Hfl, and this sug- 
gests that the NP OPE component for Q 2 < 
A 2 should be replaced by an expansion analo- 
gous to Eq.(8), but in powers of Q 2 / 'A 2 . Con- 
tinuity of cjp(t) and its first three deriva- 
tives at t — 1, and equivalently finiteness 
of V^p^{Q 2 ) and its first three derivatives 
d/dlnQ at Q 2 — A 2 is underwritten by the 
relation of Eq.(19), and by three additional 
more complicated relations involving the ^0,1 
and Bo 1 residues ^ It is easy to show that 
the ambiguous imaginary part in T>pf arising 
from IR renormalons for Q 2 > A 2 , and UV 
renormalons for Q 2 < A 2 , can be written di- 
rectly in terms of lo^ and u^. For Q 2 > A 2 
one has 



2tt A 2 



A 2 



Im[P^(Q 2 )]=±^^(^). (2:,) 

and for Q 2 < A 2 , 

,(i), n 2,_ , 2tt A 2 uv / A 2 



Im[V^{Q 2 )] = ± T ^ V [gs ) ■ (26) 

Continuity at Q 2 = A 2 then follows from con- 
tinuity of u>T>(t) at t = 1. In principle the real 
part of the OPE condensates are independent 
of the imaginary, but continuity and finite- 
ness involve the set of relations between Ao,i 
and B 1 that we have just noted. Continuity 
naturally follows if we write D^ p (Q 2 ) in the 
form 

A 2 /Q 2 



K± 



2ni 



dt uj v {t) + t 



dLUT,(t)\ 

dt ) ■ 

(27) 



Here k is an overall real non-perturbative 
constant. If the PT component is PV reg- 
ulated then one averages over the ± possi- 
bilities for contour routing, combining with 
Vp 1 ^ one can then write down the overall re- 
sult for P( L '(Q 2 ). The Q 2 evolution is fixed 
by the non-perturbative constant k and by 
A. Both the PT and NP components freeze 
smoothly to zero. 

It is a pleasure to thank the organis- 
ers of ICHEP'06 for arranging such a stim- 
ulating conference. Thanks also to Irinel 
Caprini, Jan Fischer, Andrei Kataev and 
Dimitri Shirkov for some useful and enter- 
taining discussions. 
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